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Abstract 


The OUCS2 upwind compact scheme for calculation of first derivatives in the Euler and Navier- 
Stokes equations is the focus of the present paper. Derived and analyzed by Sengupta et al. and 
primarily meant for incompressible flow problems, the scheme has recently been applied to com- 
pressible flow with shocks. To handle shocks effectively, a combination of second and fourth der- 
ivative artificial dissipation terms were used by Sengupta et al., instead of the inbuilt sixth deri- 
vative dissipation of the OUCS2. For unsteady flow problems with shocks and vortices, use of hi- 
gh order numerical dissipation is desirable for proper resolution of small scale structures. In a pr- 
esent paper we show that when used with flux limiters, one can use the inbuilt high order dissipa- 
tion in smooth regions and still bring down the scheme to small stencil formulae with low order d- 
issipation near shocks. The basic solver is built upon the AUSM™ algorithm. OUCS2 comes in w- 
hile computing the left and right states of the primitive variables at the cell faces, requiring little 
effort to modify an existing AUSM™ based solver to a high resolution version with the help of O- 
UCS2. We demonstrate the procedure for a number of problems in one and two dimensions solvi- 
ng the Euler equations. 
Keywords: AUSM*, Compressible flow, Euler equations, Flux limiter, High Mach number 


Introduction 

Numerical simulation of Euler or Navier-Stokes equa- 
tions using finite difference or finite volume based meth- 
ods requires discretization of the first and second deriva- 
tive terms in terms of values of the variables at discrete 
grid locations. Such formulae may be explicit - linking 
the derivative at a grid location to the variable value 
in a chosen neighbourhood - or they may be in the 
form of an implicit relation which includes neighbour- 
hood derivatives in a single formula. The latter class 
of schemes are known as compact schemes. Utilizing 
compact schemes for fluid dynamics applications can 
be linked to Anderi I. Tolstykh, who suggested their 
use back in 1972 (Tolstykh, 1994a). The various com- 
pact schemes derived by Tolstykh and his group be- 


long to the CUD (Compact Upwind Differencing) fam- 
ily of schemes. Developments on the third and fifth or- 
der CUD schemes may be found in (Tolstykh, 1994a,b). 
Their application to compressible flow problems are 
mentioned in (Tolstykh, 1994a, 1973; Ravichandran, 
1994, 1997; Tolstykh and Lipavskii, 1998). Application 
of compact schemes to compressible solvers demands 
special care due to the presence of sharp gradients in 
the form of shock waves. One needs to apply low or- 
der numerical dissipation locally at a shock for stable 
computation. For compressible applications, the Euler 
or Navier-Stokes equations are usually discretized in a 
conservative form involving fluxes at the cell faces. One 
may, therefore, derive the compact scheme in a form 
suitable for the flux calculation. This approach has been 
followed in (Tolstykh, 1994a; Pirozzoli, 2002). Alterna- 
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tively, the compact derivative scheme may be applied 
straightaway in a nonconservative manner, as demon- 
strated in (Sengupta et al., 2013). Though the govern- 
ing equations are in conservative form, in the latter ap- 
proach, dependence on a pressure based switch leads to 
non-cancellation of variable values in the interior, and 
boundary stencils contribute to the same effect near the 
boundaries. 


It is the demand for schemes with high spectral ac- 
curacy in smooth regions and highly dissipative schemes 
at the shock that has led to the development of a num- 
ber of methods addressing this issue. Rizzetta and Mor- 
gan, (2008) recommends use of central compact schemes 
which are nondissipative, together with compact filter- 
ing schemes. The basic solver is based on an implicit 
approximately factored finite difference scheme, and in 
the words of the authors - ‘nonlinear artificial dissipa- 
tion terms are also appended to the implicit operators 
to augment stability.’ This does not allow estimation 
of the capacity of the filtering technique Rizzetta and 
Morgan (2008) alone in the presence of strong shocks. 
An explicit solver would be ideal for the evaluation of 
the compact derivative and filtering schemes. This is 
followed by Meinke et al. (Meinke et al., 2002), who 
employ an explicit AUSM based solver together with 
compact derivative and filtering schemes. This, how- 
ever, has been tested only for large eddy simulation of 
planar turbulent jets. Kawai et al. (Kawai et al., 2010) 
use a combination of compact derivative scheme and a 
filtering scheme as above, together with a localized ar- 
tificial diffusivity (LAD) which activates near shocks. 
This method has primarily been tested for large eddy 
simulation of compressible turbulent flows. Currently 
the more popular methods for application of compact 
schemes to compressible shock dominated flows involve 
the ENO or WENO schemes. Here compact deriva- 
tives enter as part of a hybridization strategy (Pirozzoli, 
2002). In these methods, the compact scheme operates 
in the smoother regions, leaving the task of shock cap- 
turing to the ENO/WENO scheme. The general idea 
in a WENO scheme is to compute the numerical flux 
from a combination of a number of candidate stencils 
(Jiang and Shu, 1996), which renders implementation 
of the method complex. Complexity is only enhanced 
in a hybridized form with a compact scheme. 


The OUCS2 was derived and analyzed in the wavenum- 


ber space in (Sengupta et al., 2003). Being a fifth-order 
upwind compact scheme, it possesses high accuracy in 
the wavenumber space and low numerical dissipation at 
lower wavenumbers. As mentioned before, when applied 
with a combination of fourth and second derivative nu- 
merical dissipation, the scheme partly loses this attrac- 
tive property of low dissipation in the lower wavenum- 
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ber range. In order to retain its inbuilt sixth derivative 
driven numerical dissipation, we felt encouraged to ap- 
ply it in combination with a flux limiter, particularly 
after noticing the success achieved by application of lim- 
iters with high resolution explicit conservative schemes 
(Hahn and Drikakis, 2009; Halder et al., 2013; Dora 
et al., 2014; Kundu and De, 2014). The last three refer- 
ences have dealt with shock and vortex dominated flows 
simulated by an explicit AUSM™ based solver with ex- 
plicit high resolution schemes. In the following sections, 
we show that a similar strategy is also possible for the 
OUCS2 scheme - written as a cell-face interpolation 
scheme - and this yields accurate results. 
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Fig. 1. Real and imaginary components of the modified 


wavenumber of the OUCS2. For comparison, the modified 
wavenumber components of a third-order explicit upwind 
scheme (Eq. 7) is also shown. 


Numerical method 


The Euler equations for compressible flow may be writ- 
ten concisely as 


we 
at 


Of Og _ 
fest ay 


The variables can be expanded as: 


0. (1) 


p pu pv 
2 
_ | pu -_ pu’ +p 2 puv 
q= pv = puv 19 pu-+p 
pE (pE + p)u (pE + p)v 
Energy is given by pE = ®> + 5p (u? + v*). u,v are 


the velocity components in the x- and y-directions, re- 
spectively. p denotes the gas density and p denotes the 
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pressure. To comprehend the application of the com- 
pact scheme in the place of MUSCL-type cell-face inter- 
polation formula, let us have a look at how the cell-face 
flux is written in AUSM* (Liou, 1996): 


litiy = 441, 5[0.5m 442 5 (6, + Bp) — 0.5% 42 
(fp — @,)}. 
(2) 


In the above, we have omitted the pressure term 
that appears if the the above equation is for the x- 
momentum equation. Y,,1 ; = |m,41,;| in the original 
version of the AUSM?. We use a slightly modified for- 
mula which retains non-zero numerical dissipation as 
Mach number goes to zero (De and Murugan, 2011). 
The quantities 6; and ®p represent the left and right 
states composed of the left and right states of the prim- 
itive variables, as suggested by Liou (Liou, 1996). The 
OUCS2 compact scheme has been applied as an inter- 
polation formula to compute these left and right states 
of the primitive variables, as explained below. 

The OUCS2 compact upwind scheme for calculation 
of a first derivative is written as, 


/ / / 
biu;_4 + bow; + b3Uj44 = (a 2U;—2 +a 1Uj-1 





+ agu; + Q1Uj+1 + a2u;+42)/h. 
(3) 


a a 








In Eq. 3, 61 = 3 Tz 63 = 2 + 731 Oe2 = +3 + 7) 
az, = +e + 5,4 = —{, with bo = 36. a is an 
upwinding parameter. The recommended value of a is 
—40 (Sengupta et al., 2003). The real and imaginary 
components of the modified wavenumber of the OUCS2 
are displayed in Fig. 1. In the form of an interpolation 
scheme for the cell-faces, the OUCS2 scheme thus be- 


comes: 














byu;_2 + bouj42 + b3uj;43 = €_1Uj-1 + €0U; 


+ €1Uj+41 + €2Uj;+2- 








The new coefficients are e_; = be 751 CO = oe 
a — 29b2 1 a —b2 4 ao 
37 C1 = 3g + g and €2 = 36 + 7p. 


The following boundary closure is used for 
non-periodic boundaries: 


3u; + 1l2uj;41 + Uj+e2 (5) 
10 
The distortion in the modified wavenumber at the 
near-boundary points introduced by this one-sided clo- 
sure does not affect the solution as we use explicit up- 
wind schemes of third and fifth order at these points. 
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Once the cell-face values of the primitive variables 
are available from the above equation, a TVD limiter 
is applied, as is done in any high resolution solver that 
uses an explicit cell-face interpolation scheme (Kundu 
and De, 2014). The form of OUCS2 as given in Eq. 4 
makes it possible to apply the OUCS2 in a fully con- 
servative form. The CUD-3 upwind scheme described 
previously had been tested in both conservative and 
non-conservative forms (Tolstykh, 1994a). The author 
suggested avoiding the non-conservative form for com- 
puting discontinuous solutions; the conservative option 
was viewed as a ‘good choice.’ The conservative form 
of the OUCS2, suggested above, has the same advan- 
tage over its non-conservative version Eq. 3 that was 
used in (Sengupta et al., 2013). There is no necessity 
for performing any ‘symmetrization’, which was done 
before the scheme was used with numerical dissipation 
given by a combination of fourth and second derivatives 
(Sengupta et al., 2013). Very satisfactory results have 
been obtained using the OUCS2 in the form given by 
Eq. 4, described in the next section. 
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Fig. 2. Solution of the wave equation: (a) OUCS2 with 
limiter at t = 2, (b) Solution given by the a seventh-order 
WENO scheme coupled with a fifth-order compact scheme at 
t = 2. (c) and (d) are computed by the OUCS2 with limiter 
at t= 5 and 15, respectively. Number of cells = 100; CFL = 
0.5; solid lines: exact solution. 


Results and discussion 


Having outlined the method of application of the OUCS2 
scheme in an AUSMt-based Euler solver, we now focus 
on a number of standard problems in one and two di- 
mensions to assess the performance of the scheme. In 
doing so, we must remember that such a scheme of high 
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spectral accuracy is not needed near shocks. All we ex- 
pect from the solver is to be able to detect the shock 
and switch to suitable lower order formulae to avoid 
oscillatory results. Thus in the test problems we note 
whether the shock is captured without significant ‘wig- 
gles’ nearby and whether away from it the smoother 
variations in the variables are accurately represented by 
the compact scheme. We have used a four-stage Runge- 
Kutta scheme given in for time integration in all the 
problems except the first, where to facilitate compari- 
son of results we have used a different four-stage Runge- 
Kutta scheme. 





















































Fig. 3. Interaction of shock and density wave. Left frame: 
full domain, right frame: zoomed view of the fluctuating den- 
sity. Results correspond to t = 1.8. grid size: 400 cells, CFL 
= 0.1. 


Problems in one dimension 


We will first solve a linear one dimensional wave prop- 
agation problem to demonstrate how the accuracy of 
the compact scheme can be used to our advantage if 
a proper limiter is in action. Further, two one dimen- 
sional standard test cases are computed solving the Eu- 
ler equations. 


Wave propagation in one dimension 


The governing equation is 


du 
at 


Ou 


a 6 
a (6) 
The initial condition that is specified at t = 0 prop- 
agates to the right with a speed of c = 1. We specify 
the following discontinuous initial condition: 


0. 


—x sin (3727) 5 -l<a< —; 
u(x + 0.5) = ¢ |sin (272) |, |x| < 3 
2n —1—sin(3rx) /6. ¢<2<1 


Periodic boundary conditions are applied at both 
the ends of the domain. This problem was solved in (Shu 
and Osher, 1989). In Fig. 2. we compare the result of 
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the OUCS2 at t = 2 with (Pirozzoli, 2002), both com- 
puted using 100 cells and the same Runge-Kutta time 
stepping. At t = 2 the initial wave travels a distance 
of one wavelength, and appears at the same initial po- 
sition. In (Pirozzoli, 2002) this problem was computed 
using a seventh-order WENO scheme coupled to a fifth- 
order compact upwind scheme. It has been shown in 
(Sengupta et al., 2004) that for such problems it is suf- 
ficient to compute the first derivative using a ‘strictly 
one-sided’ explicit third-order upwind scheme near dis- 
continuities. The discontinuities were detected by the 
magnitude of the second derivative. When | ex > 50, 
the one-sided formula was used, yielding results that 
were superior to those obtained by the more complex 
WENO-compact hybrid scheme referred to above. How- 
ever, the comparison was made using 200 cells, citing 
greater difficulty in computing with less dissipation as- 
sociated with more number of cells. We choose to com- 
pare with the results in (Pirozzoli, 2002) using the same 
100 cells in the domain. Though not shown here, with 
200 cells the OUCS2 produces even better results in 
our mode of application, contrary to the reverse trend 
noticed in (Sengupta et al., 2004). For wave propaga- 
tion phenomena in one dimension with discontinuous 
initial data, it is sufficient to locate the presence of the 
discontinuities and apply a low-order dissipation. We 
detect the discontinuities by using Harten’s TVD con- 
dition in the limiter, and apply a second derivative as 
artificial dissipation in a conservative manner. Thus we 
are not using any problem-dependent threshold value, 
as was done in (Sengupta et al., 2004), and the TVD 
condition violation proves sufficient to detect the dis- 
continuities. We also do not switch to a fully explicit 
one-sided form of the first derivative, and the compact 
scheme is retained near discontinuities - only extra low- 
order numerical dissipation is added. The Fig. 2. shows 
that at t = 2 (see frames (a) and (b)) the OUCS2 pro- 
duces better results compared to the compact-WENO 
scheme. 
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Fig. 4. Interaction of blast waves in one dimension: result 
at t = 0.038. Left frame: full domain, right frame: zoomed 
view. Grid size: 800 cells, CFL = 0.1. 
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For this kind of problems it is possible to obtain a 
good match with the exact solution momentarily while 
the solver actually is on the way to produce oscillatory 
diverging results at greater time instants. To show that 
this is not the case with the present form of application 
of the OUCS2, the frames (c) and (d) are presented, 
which show the results at t = 5 and 15. We notice that 
the time iterations continue in a stable manner without 
any new extrema. The effect of inbuilt artificial dissi- 
pation of the overall scheme is reflected in reduction of 
the amplitude with time. The results for this test case 
proves that in one-dimensional wave propagation prob- 
lems, use of compact schemes in a conservative form, as 
in the present case, can be very accurate in combination 
with a TVD limiter. 











Fig. 5. Shock diffraction at a corner: numerical schlieren at 
t = 1. OUCS2 scheme on 1000 x 1000 grid. Computed with 
a constant At = 10~+. 40 contour levels have been plotted 
from 1 to 40 at an increment of 1. 


Shock-density wave interaction 


This problem is taken from (Shu and Osher, 1989). A 
Mach 3 shock interacts with a density wave in a domain 
—5 <a <5. The initial condition is 


iomn)={ 


400 cells are taken in the domain. We solve the Eu- 
ler equations in one dimension and compare results of 
OUCS2 with the following standard third-order explicit 
upwind scheme: 


(3.857143, 2.629369, 10.333333) , 
(1+ 0.2sin (52) ,0, 1), 


a<—A4 
z>—A 


_ QUi41 + DU; — Uj-1 


; (7) 





Uy 


NIR 
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In Fig. 3. the result at t = 1.8 is shown. The refer- 
ence solution is computed on a grid of 4000 cells using 
a ninth-order explicit scheme. We note that the den- 
sity fluctuation is better captured by the OUCS2 com- 
pared to the third-order explicit scheme. This shows 
that the limiter is active only in shock-capturing, and 
the OUCS2 is resolving the smoother fluctuations. 














Fig. 6. Shock diffraction at a corner: numerical schlieren at 
t = 1. third-order explicit scheme (Eq. 7) on 1000 x 1000 grid. 
Computed with a constant At = 10-4. Same contour levels 
as in Fig. 5. 


Interacting blast waves 


For this problem the initial condition is (Shu and Osher, 
1989) 


(1,0,1000), 0<2<01 
(p,u,p) = § (1,0,0.01), O1<2<0.9 
(1,0,100), 0.9<a<1 


We solve the Euler equations in one dimension us- 
ing 800 cells and view the result in Fig. 4. Solid wall 
boundary condition is applied at the ends of the do- 
main. This problem was chosen to highlight the fact 
that the present compact scheme-limiter combination 
is able to handle a pressure jump by a factor of 10° 
in the initial conditions. Density variation is generally 
sharper in the OUCS2 result which is visible in the 
sharp density rise at « = 0.6. If we compare with the 
results in figure 13 of (Shu and Osher, 1989), we note 
that the OUCS2 result is comparable to the result of 
the ENO-RF-S-3 computed on the same grid. 
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Problems in two dimensions 


We now focus on three two dimensional problems which 
involve both shock waves and vortices. There are prob- 
lems where high wavenumber excitation is nearly ab- 
sent. In such cases a high accuracy compact scheme and 
a lower order explicit scheme produce comparable re- 
sults. We choose our test cases carefully so that the use 
of a compact upwind scheme makes a difference when 
compared with a third-order explicit formula. 


Shock diffraction at a corner 


This problem is taken from (Artebrant and Schroll, 
2006). The domain of interest is a unit square. Initially, 
the whole of the interior has u= v =0 and p=p=1. 
At the inflow, which is the top half of the left verti- 
cal boundary, p = 3 and flow angle is 7/10. Solid wall 
boundary conditions are set at the lower half of the 
left boundary and the bottom wall. Artificial absorbing 
boundary condition is set at the top and right walls. 
This is similar to the Schardin’s problem (Halder et al., 
2013), but the present domain is small and set-up of 
the problem much easier. 





Fig. 7. The four-shocks problem: numerical schlieren at ¢ = 
0.8. (a) OUCS2 on 1000 x 1000 grid, (b) third-order scheme 
on 1000 x 1000 grid, (c) OUCS2 on 2400 x 2400 grid, (d) 
third-order scheme on 2400 x 2400 grid. CFL = 0.5 for all. 40 
contour levels with an increment of 1. 


Note that the initial condition is not fully correct in 
(Artebrant and Schroll, 2006). The shock wave that is 
part of the initial condition diffracts at the corner 
(a =0, y = 0.5 - the point separating the inflow and 
the wall boundary condition at the left vertical bound- 
ary). 

The diffracted shock is later reflected back into the 
domain at the bottom wall. Meanwhile, a vortex grows 
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behind it near the inflow. When the reflected shock 
moves upward, it interacts with this primary vortex. 
Small-scale vortices are also noticed behind the pri- 
mary vortex. These are generated due to the Kelvin- 
Helmholtz instability. In Fig. 5., we notice these small- 
scale vortices moving into the primary core. The re- 
flected shock interacts with the primary vortex. It can- 
not pass through the vortex with its shape intact, and 
breaks into accelerated and decelerated components - 
visible in Fig. 5. 

In Fig. 6., the solution produced by the third-order 
scheme (Eq. 7) is shown. Comparing, we note that the 
OUCS2 result has more small-scale trailing vortices. 
This is indicative of the high wavenumber accuracy of 
the OUCS2. Similar vortices are also generated behind 
a shock-tube generated vortex ring (Murugan et al., 
2012). Using a second-order accurate scheme, these could 
not be detected as early as they could be seen in exper- 
iments. OUCS2 or a similar scheme is likely to detect 
these small-scale vortices better. The background dis- 
turbance in the OUCS2 result occurs for other high 
resolution schemes also. This has been tested by us us- 
ing another high resolution scheme. This is likely to be 
the result of artificial inflow boundary condition applied 
here. This disturbance remains undetected by third or- 
der accurate schemes. 











Fig. 8. The four-shocks problem: numerical schlieren at 
t = 0.8. OUCS2 on 2400 x 2400 grid. High wavenumber accu- 
racy of the OUCS2 generates the smaller-scale structures that 
would not be visible by lower accuracy schemes. 50 contour 
levels are drawn between 1 to 99 at an increment of 2. 


The four-shocks problem 


In this problem, the initial condition contains four shocks 
- dividing the unit square domain into four quadrants. 
The variables are set as 
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(1.5, 0,0, 1.5), Ql 
(0.53225806, 1.20604538, 0, 0.3) , Q2 
(p, u,v, p) = & (0.13799283, 1.20604538, 1.20604538, 
2.90322581 x 1072), Q3 
(0.53225806, 0, 1.20604538, 0.30) , Q4 


where Q1 : [0.75, 1] x [0.75, 1]; Q2 : [0, 0.75] x [0.75, 1]; 
Q3 : [0, 0.75] x [0,0.75] and @Q4 : [0.75, 1] x [0, 0.75]. 

The full description is provided in (Schulz-Rinne 
et al., 1993). We keep the point of contact of the four 
shocks at (2 = 0.75, y = 0.75), instead of at the cen- 
tre of the domain. This was also done in (Serna, 2006) 
for better visualization of the results. We did not type 
in the values given above, rather (Schulz-Rinne et al., 
1993) was consulted for analytical expressions when- 
ever available. With time, four slip lines develop. In 
Fig. 7. the results of the OUCS2 and the third-order 
explicit scheme, given by Eq. 7, are plotted for two dif- 
ferent grid sizes - 1000 x 1000 and 2400 x 2400. Both 
the OUCS2 and the third-order scheme produce the 
key structures similar to the ones in published litera- 
ture (Schulz-Rinne et al., 1993; Serna, 2006; Cada and 
Torrilhon, 2009). The last reference solved this prob- 
lem on a 1000 x 1000 grid using a second-order accu- 
rate TVD-MUSCL scheme with various limiters. For a 
certain limiter (designated LimO3), the roll-up of the 
slip lines was most prominent. 
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Fig. 9. The four-contacts problem: numerical schlieren at 
t = 0.8. (a) OUCS2 and (b) third-order scheme on 2400 x 2400 
grid. CFL = 0.5. 50 contour levels are drawn between 1 to 99 
at an increment of 2. 


In our case, this roll-up is much stronger for the 
OUCS2 scheme on the same grid. The third-order scheme 
also produces roll-up structures, but to a less extent. On 
the 2400 x 2400 grid, the OUCS2 creates small scale 
details which are resolved only due to its high spectral 
accuracy for a larger wavenumber range. This is shown 
more clearly in Fig. 8., which shows a zoomed view of 
the central region. The finer details visible in this figure 
confirm the suitability of the OUCS2 for shock-vortex 
interaction cases where small-scale structures are gen- 
erated. 
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The four-contacts problem 


The initial condition for the four-contacts problem is 
given by (Schulz-Rinne et al., 1993) 


1.0, 0.75, —0.50, 1.0) , Ql 
2.0, 0.75, 0.50, 1.0) , Q2 
1.0;—0.75,0.50;1.0 x 10"), 93 


( 

_ JC 

(p, U, U, P) = ( 
(3.0, —0.75, —0.50, 1.0) , Q4 


where Q1 : [0.5,1] x [0.5,1]; Q2 : [0,0.5] x [0.5, 1]; 
Q3 : [0,0.5] x [0, 0.5] and Q4 : (0.5, 1] x [0, 0.5]. 











Fig. 10. The four-contacts problem: numerical schlieren at 
t = 0.8. (a) OUCS2 on 2400 x 2400 grid: zoomed view to show 
the small-scale structures. Same contour levels as in Fig. 9. 


The solution is a clockwise turning vortex with four 
slip lines spiralling around the centre (Schulz-Rinne et al., 
1993). For this problem, the conditions listed above can 
be typed in, and no analytical expression is involved for 
accurate representation. Here also the slip lines gener- 
ate roll-up vortices. Results are plotted in the form of 
numerical schlieren at t = 0.8 in Fig. 9. From our sim- 
ulation results and the ones in (Cada and Torrilhon, 
2009), we note that the slip lines which extend to the 
first and the third quadrants are the ones which develop 
the smaller scale vortices. In fact, these slip lines do not 
show any such vortices prominently in (Cada and Tor- 
rilhon, 2009) for the second-order TVD-MUSCL scheme 
with any of the limiters used on a 1000 x 1000 grid. Only 
the slip lines extending to the second and the fourth 
quadrants develop a couple of roll-up vortices. In our 
case a similar trend is observed. In Fig. 9., the slip lines 
in the second and the fourth quadrants display these 
roll-up structures all along their length for both OUCS2 
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and the third-order scheme. In the other two quadrants, 
only OUCS2 is able to excite these structures along the 
whole length of the slip lines. The third-order scheme 
does not show these near the outer ends. Fig. 10. is a 
zoomed view of the central region for OUCS2. The inte- 
rior of the rolled-up structures now show even smaller- 
scale features. The third-order scheme is unable to pro- 
duce them, and this is where a fifth-order upwind com- 
pact scheme like the OUCS2 can be most effectively 
used. 


Concluding remarks 


We have presented a method of application of the OUCS2 
upwind compact scheme for compressible flow prob- 
lems with shocks. The OUCS2 is a fifth-order upwind 
compact scheme with artificial dissipation provided by 
a sixth-derivative. Our approach enables applying the 
OUCS2 scheme to its full strength, and not by applica- 
tion of lower-order (fourth) dissipation explicitly, as has 
been done in recent literature. We apply the OUCS2 in 
a conservative manner, instead of losing discrete con- 
servation due to its use in the derivative form cou- 
pled with switches to identify shocks. No symmetriza- 
tion is performed, and the OUCS2 is integrated into an 
AUSM*-based Euler solver equipped with TVD lim- 
iting. A number of test cases in one and two dimen- 
sions indicate that accurate results can be obtained by 
the present approach. For one-dimensional linear wave 
propagation with discontinuous initial condition, the 
OUCS2 combines well with a TVD limiting strategy 
to capture the solution more accurately in comparison 
to a seventh-order WENO scheme hybridized with a 
fifth-order compact scheme. Two other standard one- 
dimensional problems of gas dynamics show the accu- 
racy of the OUCS2 in smooth regions and its ability to 
combine with the TVD limiting procedure to predict 
evolution of sharp jumps in the initial condition. The 
two-dimensional problems have small-scale vortices and 
shocks/slip lines developing into an unsteady flow-field. 
The OUCS2 has been shown to be effective in capturing 
these small-scale structures much better than a stan- 
dard third-order explicit scheme. Such structures are 
part of many shock-vortex and shock-turbulence inter- 
action problems, and this is where we plan to apply the 
OUCS2 and other high spectral accuracy schemes in 
future. 
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